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SUMMARY 


NAE  and  de  Havilland  Aircraft  of  Canada  have  designed  and 
tested  in  the  NAE  5  ft  X  5  ft  wind  tunnel  a  21%  thick  supercritical  airfoil 
capable  of  sustaining  long  runs  of  laminar  flow  on  both  surfaces.  The 
measured  drag  levels  are  superior  to  those  of  any  model  previously  tested  in 
this  facility  and  are  comparable  to  classical  NACA  and  modem  NASA  NLF 
airfoils. 


RESUME 


NAE  et  de  Havilland  Aircraft  of  Canada  ont  congu  et  teste  dans  la 
soufflerie  aerodynamique  NAE  5  pi  X  5  pi  un  profil  supercritique  d’une 
epaisseur  de  21%  capable  de  maintenir  de  longs  trajets  d’ecoulement 
laminaire  sur  ses  deux  surfaces.  Les  niveaux  de  trainee  mesures  sont 
superieurs  a  ceux  de  tous  les  modeles  testes  precedemment  dans  ces  instal¬ 
lations  et  sont  comparables  aux  profils  NACA  classiques  et  NASA  NLF 
modemes. 
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SYMBOLS 

Definition 
Mach  number 
lift  coefficient 

lift  coefficient  at  which  separation  first  starts  to  occur 

pitching  moment  about  quarter  chord,  negative  nose  down 

drag  measured  by  the  wake  rake  and  averaged  from  probes  1  and  3  (see  Fig.  10) 

drag  measured  from  probe  2 

drag  measured  from  probe  3 

Reynolds  number  based  on  a  10  in.  chord 

angle  of  attack  (corrected  for  wall  interference) 

drag  rise  Mach  number 

taken  from  balance  measurements 

free  stream  (free  stream  Mach  number  is  corrected  for  wall  interference) 


max 


maximum  value 
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ANALYSIS  OF  EXPERIMENTAL  DATA  FOR  A  21%  THICK  NATURAL 
LAMINAR  FLOW  AIRFOIL,  NAE  68-060-21: 1 


1.0  INTRODUCTION 

As  part  of  an  ongoing  study  of  supercritical  wing  sections  carried  out  jointly  by  NAE  and 
de  Havilland*,  we  present  here  some  results  obtained  from  wind  tunnel  tests  on  a  21%  thick  airfoil. 
Other  airfoils  in  the  series  are  16%  thick  (see  Ref.  1)  and  shortly  will  include  13%  and  10%  thick 
sections.  The  main  objective  of  the  present  design  is  to  investigate  the  possibility  of  achieving  long 
runs  of  laminar  flow  on  both  the  upper  and  lower  surfaces.  It  was  thought  this  would  be  difficult  in 
the  NAE  5  ft  X  5  ft  blowdown  wind  tunnel  since  it  is  not  a  low  turbulence  tunnel  and  the  effect  of 
the  free  stream  turbulence  on  transition  was  uncertain. 

The  study  of  Natural  Laminar  Flow  (NLF)  airfoils  has  been  in  progress  worldwide  for  many 
years  but  it  has  always  been  difficult  to  achieve  good  NLF  results  in  real  flight  due  to  imperfections  in 
the  wing  surfaces.  Recently  there  has  been  a  more  sustained  interest  in  the  NLF  concept,  particularly 
in  the  United  States,  as  modem  manufacturing  methods  and  materials  yield  much  improved  wing 
surfaces.  Unfortunately  most  of  the  modem  work  is  classified  and  so  comparisons  of  our  own  work  in 
this  field  with  ethers  is  difficult.  Only  at  low  speed  (Somers,  Refs.  2  and  3)  are  other  results  available 
and  comparisons  are  made  in  Reference  1  which  show  that  our  drag  levels,  at  supercritical  Mach 
numbers,  are  very  similar  to  drag  at  low  speed,  reported  by  Somers,  in  the  Langley  low  turbulence 
tunnel. 


The  NLF  concept  for  achieving  low  drag  is  only  one  method  amongst  many  others  such  as 
suction  and  energizing  of  boundary  layers  on  aerodynamic  surfaces,  Reference  4,  wavy  surfaces. 
Reference  5,  augmentor  airfoils,  Reference  6,  large  eddy  break-up  devices.  Reference  7,  turbulence 
manipulators.  Reference  8. 

The  computational  method  employed  in  our  design  was  the  well  known  BGK  Computer 
code!9!  modified  to  include  Green’s  boundary  layer  method.  The  non-conservative  option  was  used 
as  this  gives  a  better  correlation  with  experiment.  A  slightly  favourable  but  fairly  flat  upper  surface 
pressure  distribution  up  to  the  shock  was  sought  at  the  design  condition  (M^  =  0.68,  CL  =  0.6)  while 
on  the  lower  surface  a  quite  favourable  gradient  up  to  minimum  pressure  at  about  40%  chord  was 
sought.  A  modest  aft  loading  was  also  a  criterion  but  with  pitching  moments  not  too  high.  These 
conditions  were  met  using  the  BGK  code  with  transition  fixed  at  the  standard  7  and  15%  locations 
on  the  upper  and  lower  surfaces.  It  is  demonstrated  that  a  reasonable  agreement  is  obtained  between 
this  theoretical  pressure  distribution  and  that  obtained  from  the  experiment.  In  the  experiment  one 
must  remember  that  the  pressure  holes  themselves  cause  turbulence  and  so  the  measured  pressure 
distribution,  in  a  turbulent  chordwise  strip,  is  presumably  not  representative  of  the  NLF  pressure 
distribution  on  the  remainder  of  the  airfoil.  Thus  any  prediction  of  the  NLF  pressure  distribution 
cannot  be  verified  under  present  conditions  but  would  need  unobtrusive  measurement  techniques. 
Alternatively  two  models,  one  with  staggered  pressure  holes,  could  be  used.  The  BGK  prediction  of 
pressures  with  a  free  transition  option  developed  by  de  Havilland  shows  a  more  favourable  distribution 
on  the  upper  surface  and  a  less  favourable  one  on  the  lower  surface  compared  to  7  and  15%  transition 
points.  It  also  shows  more  aft  loading  but  the  same  shock  strength. 

Our  studies  also  show  that  as  the  free  stream  Mach  number  increases  through  supercritical, 
at  the  design  lift  say,  the  upper  surface  pressure  gradients  become  more  favourable  as  the  shock  gets 
stronger.  Thus  it  seems  that  more  and  more  laminar  flow  is  possible  on  the  upper  surface  yielding 
lower  drags  until  eventually  the  wave  drag  dominates  the  situation  and  drag  rise  appears.  Hence  most 
of  our  drag  plots  show  fairly  extensive  buckets  covering  a  AM  of  up  to  0.04.  They  also  show  a  sharp 
drag  rise  after  the  bucket  which  may  be  due  to  laminar  flow  separation. 
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It  is  also  shown  that  fully  turbulent  drag  levels  (one  of  the  wake  rake  probes  was  directly 
in  line  with  the  pressure  tappings)  are  of  a  reasonable  size  and  are  quite  comparable  with  other  super¬ 
critical  airfoils. 

Other  experimental  data  on  thick  supercritical  airfoils  is  rather  scarce.  Van  Egmond  and 
Rozendalt10!  report  on  an  18%  thick  airfoil  (NLR  7501)  tested  in  the  NLR  pilot  tunnel  at  Mach 
numbers  of  about  0.73  but  at  low  Reynolds  number  (about  2  million).  In  terms  of  drag,  with  free 
transition,  they  found  a  narrow  drag  bucket  at  M^,  =  0.765  (uncorrected)  with  CL  ^  0.4  which 
is  the  ‘tunnel  design  condition’.  The  drag  here  was  about  100  counts  rapidly  rising  to  120  counts 
for  AMqo  —  ±0.01.  The  authors  assessed  the  narrow  bucket  as  being  due  to  flow  separation  for  a 
slight  variation  from  the  design  condition.  The  design  of  this  airfoil  was  not  favourable  to  long  runs 
of  laminar  flow  with  pressure  peaks  on  the  upper  and  lower  surfaces  being  at  about  7%  and  6% 
respectively. 

Blackwell  (Ref.  11)  presented  data  on  supercritical  airfoils  of  10, 16  and  21%  maximum 
thickness.  He  was  particularly  looking  at  Reynolds  number  scale  effects.  With  regard  to  the  21% 
airfoil,  suitable  for  span  loader  type  aircraft,  data  was  presented  for  Reynolds  numbers  of  7, 11  and 
22  million.  For  a  Mach  number  of  0.68  and  CL  about  0.65  it  was  shown  that  the  drag  decreased  from 
about  0.0183  to  0.0163  as  Reynolds  number  increased  from  7  to  22  million.  In  this  case  transition 
was  probably  at  about  10%  on  the  upper  surface  and  40%  on  the  lower  surface  judging  from  the 
pressure  distribution.  Quite  strong  shocks  were  present  on  the  upper  surface  and  the  flow  went  super¬ 
critical  on  the  lower  surface.  Blackwell  states  this  was  not  an  optimum  design  and  later  Lockheed 
developed  a  20%  thick  optimized  airfoil  which  is  classified  l12' . 

In  our  21%  airfoil  case  it  will  be  seen  that,  contrary  to  the  above,  drag  increases  over  the 
Reynolds  number  range  covered,  that  is  from  7  X  106  to  17  X  106.  This  is  presumably  due  to  the  fact 
that  we  have  a  fairly  constant  transition  Reynolds  number  so  that  transition  will  be  further  forward 
on  the  foil  as  Rc  is  increased.  Thus  there  is  less  laminar  flow  and  hence  higher  drag  as  Rc  increases. 

The  design  conditions  for  the  present  foil  were  Mw  =  0.68  with  CL  =  0.6.  Three  Reynolds 
numbers  (based  on  a  10  inch  chord)  of  6.8, 12.8  and  16.7  were  used  in  the  experimental  tests  which 
were  carried  out  in  the  NAE  5  ft  X  5  ft  blowdown  wind  tunnel.  A  data  report  (Ref.  13)  describes  the 
operating  conditions,  the  equipment  and  the  measurements  taken  in  this  test.  The  analysis  of  the 
measurements  is  the  subject  of  this  report. 

2.0  LIFT  CHARACTERISTICS 
2.1  A  Discussion  of  CL-a 

In  all  our  discussions  of  lift  and  moments  we  decided  to  use  balance  measurements  as  these 
better  represent  laminar  flow  conditions  rather  than  pressure  integrated  values.  In  actual  fact  the  two 
readings  did  not  give  significantly  different  values. 

Two  types  of  graph  are  presented  in  this  discussion.  In  the  first  set,  Figures  1  and  2, 
Reynolds  number  Rc  (based  on  chord)  is  held  constant  whilst  the  Mach  number  is  changed.  In  the 
second  set,  Figures  3  and  4,  the  converse  is  true. 

The  pitching  moment  shown  in  the  above  figures  will  be  dealt  with  later. 

The  constant  Reynolds  number  graphs  of  Figures  1  and  2  show  that  increasing  the  Mach 
number  has  the  effect  of  increasing  dCL/da  in  the  linear  part  of  these  curves.  It  changes  gradually 
from  a  low  value  of  d CL/da  =  0.116  at  M^,  =  0.3  to  dCL/da  =  0.2  at  M^  =  0.7.  Figures  1  and  2  also 
show  that  in  order  to  obtain  the  design  Clb  =  0.6  at  constant  Rc  =  12.7  X  106,  in  the  Mach  number 
range  0.3  <  M^,  <  0.7,  the  angle  of  attack  must  correspondingly  reduce  from  1.718  to  0.5°.  The  zero 
incidence  lift  coefficient  Clq  _q  also  increases  with  increasing  Mach  number  from  a  low  of  about 
Clb  =  0.4,  at  M^,  =  0.3  to  a  peak  of  Clb  =  0.513  at  M^  =  0.68. 
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The  three  plots  in  Figures  3  to  5  show  the  effect  of  changing  the  Reynolds  number  whilst 
the  Mach  number  is  held  constant.  It  seems  that  on  the  linear  portion  of  the  Clb  versus  a  curve 
the  Reynolds  number  change  has  small  effect,  except  at  =  0.7  (Fig.  5)  where  the  lower  Rc  case 
produces  significantly  more  lift.  It  can  be  seen  from  these  figures  that,  as  expected.  Cl  is  quite  linear 
up  to  stall  onset.  Beyond  this  point  the  stall  is  well  behaved  and  does  not  produce  any  rapid  loss  of 
lift. 


The  CLmax  obtained  is  very  much  a  function  of  the  Mach  number.  The  highest  CLmax  value 
obtained  was  about  1.7  at  M^  =  0.3,  Rc  =  6.84  X  106  at  a  =  13.75°  shown  in  Figure  3.  At  the  design 
Mach  number  M^,  =  0.68,  and  Rc  =  12.87  X  106  the  CLmax  obtained  was  about  0.95  at  a  =  3.4°  a*, 
seen  in  Figure  2.  By  relaxing  the  Reynolds  number  down  to  Rc  =  6.83  X  106,  Cl  reaches  a  maximum 
value  of  about  1.0  at  a  =  7.4°,  see  Figure  6. 

Finally  on  Figure  7  we  show  the  variation  of  bC^ldoi  against  free  stream  Mach  number. 
These  values  are  practically  identical  for  both  Rc  =6.8  and  12.8  X  106 ;  there  is  insufficient  data  at 
Rc  =  16.7  X  106.  Also  shown  is  the  Prandtl  Glauert  curve  for  dCi/da. 


2.2  Lift  Performance 

Figures  8  and  9  show  the  lift  performance  against  Mach  number  at  low  (Rc  =  6.8  X  106) 
and  high  Reynolds  number  (Rc  =  12.8  X  106,  16.7  X  10°)  respectively.  Hie  upper  curve  in  both 
figures  corresponds  to  the  CLmax  obtained  from  the  lift-incidence  curves,  whilst  the  lower  curve 
corresponds  to  CLse  .  CLsep  is  obtained  by  a  de  Havilland  procedure  of  plotting  Cp  versus  Cl  at  a 
chord  station  x/c  =  0.96,  and  then  determining  the  point  on  the  curve  where  6Cp/90.  =  -0.4 
(Ref.  14).  Lift  corresponding  to  the  drag  rise  Mach  numbers  was  also  determined  from  appropriate 
drag  polars,  as  explained  later,  and  is  shown  in  both  figures. 

In  both  figures  it  is  apparent  that  CLmix  and  Cls?p  generally  diminish  as  Mach  number 
increases.  At  the  lower  Reynolds  number  Rc  =  6.8  X  106,  (Fig.  8)  the  rate  of  decay  of  CLmax  and 
Cl  from  a  high  of  1.69  and  1.58  respectively  at  M^  =  0.3  to  a  low  of  0.85  and  0.48  is  less  orderly 
compared  to  the  high  Reynolds  number  case  in  Figure  9.  Note  the  large  gap  between  Cl;i1v,  and  Clscp 
at  Mqo  =  0.5  in  Figure  8.  Both  Reynolds  numbers  show  a  converging  of  CLmax  and  Cj  ir»  the  region 
0.64  <  Moo  <  0.66.  At  design  Mach  number  Mw  =  0.68,  both  figures  show  values  of  CLmax  =  1.0  and 
CL  =  0.9. 

Lsep 

Very  little  data  was  available  for  Rc  =  16.7  X  106. 


3.0  A  DISCUSSION  OF  CM  -  a 

At  low  Mach  number  in  the  range  0.3  to  0.6,  Rc  =  12.7  X  106,  the  pitching  moment  Cmb 
(referenced  to  Vi  -  chord)  shows  very  small  variations  with  a,  Figure  1.  For  the  range  where  Clb  is 
less  than  CLmax  the  value  of  Cmb  lies  in  the  range  0.11  <  |Cmb  I  <  0.135.  At  the  higher  Mach  number 
range  M,*,  =  0.66  to  0.7  in  Figure  2,  there  seems  to  be  somewhat  of  a  linear  trend  of  Cmb  against  a  in 
the  range  -2  <  a0  <  1.50.  The  SCm/So;  value  is  about  -0.005,  yielding  an  aerodynamic  centre  of 
about  0.03  i.e.  28%  chord. 


At  Mach  numbers  0.3  and  0.6,  in  Figures  3  and  4,  respectively,  the  pitching  moment  is  not 
affected  significantly  by  changing  the  Reynolds  number  from  6.8  X  10°  to  12.8  X  106.  At  the  higher 
Mach  number  M^  =  0.7  (Fig.  5),  the  difference  in  Cmb  values  for  the  same  Reynolds  number  change 
is  quite  measurable.  The  values  of  Cmb  at  the  lowest  Mach  number  M,*,  =  0.3  in  Figure  3  remain  well 
behaved  within  the  range  0.1  <  |CmbI  <  0.125.  At  higher  Mach  numbers  Mm  =  0.6  and  0.7  in 
Figures  4  and  5  respectively,  Cmb  displays  a  weak  sinusoidal  relationship  with  a.  The  static  damping 
coefficient  dC\iB/da  has  values  of  -0.006  for  M^  =  0.6  (Figure  4),  and  -0.01  at  M^  =  0.7  (Figure  6) 
for  small  angles  of  attack. 


4.0  AIRFOIL  DRAG 


4.1  A  Discussion  of  Wake  Drag 

The  wake  drag  was  measured  by  the  standard  sidewall-mounted  traversing  rake  supporting 
four  pitot  probes.  The  method  is  based  on  measuring  the  momentum  defect  in  the  wake,  see 
Reference  15  for  instrumentation  and  method.  The  signal  from  probe  4  is  ignored  as  it  is  sometimes 
affected  by  the  disturbed  sidewall  boundary  layer.  Ordinarily  the  total  effective  drag  is  computed 
from  an  average  of  the  remaining  three  probes.  However,  in  this  experiment  the  second  probe 
(probe  2)  lined  up  exactly  downstream  of  the  chordwise  pressure  taps  on  the  model  surface  (see 
Fig.  10).  This  was  at  first  thought  to  be  a  bad  choice  in  that  one  probe,  with  reading  from  a  turbulent 
strip,  would  be  useless.  However  on  second  thoughts  it  was  realized  that  the  readings  could  be  useful 
for  giving  turbulent  drag  levels  more  typical  of  a  conventional  aircraft  with  imperfection  on  the  wings. 
Some  background  of  our  design  ideas  here  will  be  useful. 

The  design  of  the  21%  airfoil  was  done  in  such  a  way  as  to  capitalize  upon  long  runs  of 
laminar  flow  on  both  the  upper  and  lower  surfaces.  This  was  achieved  by  pressure  gradients  favourable 
enough  to  overcome  the  slight  imperfections  or  dust  always  present  on  a  model  surface.  To  aid  the 
latter  problem  the  model  was  cleaned  between  runs.  However,  as  mentioned,  the  pressure  taps  them¬ 
selves  gave  a  problem  since  they  produced  premature  transition  and  turbulence. 

With  the  above  in  mind  a  comparison  of  drag  from  Probe  2,  CDw2>  is  made  against  the 
average  drag  from  Probes  1  and  3  denoted  Cdw-  As  expected  the  former  is  larger  —  often  by  a  con¬ 
siderable  amount.  This  is  illustrated  in  Figures  11  to  15.  In  particular  the  ‘fully  turbulent’  drag  is 
50  counts  higher  than  the  ‘natural  laminar  flow’  drag  at  design  conditions  (M^  =  0.68,  Cl  =  0.6). 

Figures  16a  and  16b,  respectively  show  CdW2  and  Cdw  against  Mach  number  for  Reynolds 
numbers  (Rc)  6.8  X  106,  12.8  X  106  and  16.7  X  106,  at  CL  =  0.3.  CDw2  in  Figure  16a  shows  a 
continuously  increasing  trend  against  Mach  number  up  to  M^  =  0.64,  after  which  the  low  Reynolds 
number  (Rc  =  7  X  106 )  curve  begins  to  du.iinish  describing  a  bucket  with  a  minimum  value  of  about 
Cdw  =  0.0110  at  Mqq  =  0.685.  The  higher  Reynolds  number  curve  maintains  its  upward  trend, 
increasing  sharply  from  a  value  of  CpW2  =  0.013  at  M^  =  0.66  to  CpW2  =  0.017  at  M^  =  0.68  and 
then  dropping  back  to  about  CdW2  =  0.0148  at  M,*,  =  0.7.  The  reason  for  this  behaviour  is  not 
apparent.  At  Rc  =  16.7  X  106  only  three  points  were  measured,  and  they  indicate  an  increase  in  CpW2 
from  0.0108  to  0.0130  between  M^,  =  0.66  and  0.7. 

The  corresponding  drag  as  obtained  from  Cdw  for  the  above  three  Reynolds  numbers  at 
Cl  =  0.3  is  smaller  and  shows  more  distinctly  a  bucket,  see  Figure  16b.  Generally  the  curves  indicate 
a  decrease  in  drag  coefficient  up  to  M^  =  0.6.  Beyond  this  point,  there  is  a  mild  increase  up  to 
Moo  =  0.66  after  which  the  bucket  phenomenon  is  observed.  The  drag  at  the  low  Reynolds  number 
(Rc  =  6.8  X  106)  dips  as  low  as  Cqw  =  0.0057  at  M^  =  0.68.  At  Rc  =  12.8  X  106,  this  minimum 
bucket  point  has  a  value  of  about  Cdw  =  0.0065  at  Mm  =  0.685.  At  the  high  Reynolds  number 
Rc  =  16.7  X  106,  the  three  data  points  also  demonstrate  the  bucket  effect  giving  a  minimum 
CDw  =  0.0099  at  Mm  =  0.68. 

Two  approaches  are  adopted  to  evaluate  the  drag  rise  Mach  number,  MDR .  One  method  is 
based  on  determining  the  point  on  the  Drag  versus  Mach  number  graph  where  3Cd/9M  =  0.1.  The 
second  method  involves  adding  20  drag  counts  to  the  average  of  the  lower  Mach  number  Cpw  values, 
and  finding  the  point  on  the  rising  branch  of  the  curve  which  corresponds  to  this  drag  value.  The  drag 
rise  Mach  number  on  each  curve  for  the  above  two  methods  is  appropriately  shown.  The  two  methods 
give  different  values  of  MDR  as  can  be  judged  by  inspecting  any  one  curve.  For  instance,  MpR  based 
on  9Co/dM  =  0.1  for  Rc  =  6.8  X  106  in  Figure  16b  is  0.708,  and  based  on  the  20  counts  approach 
Mdr  is  0.720. 


As  expected,  the  two  drags,  (C^  and  Ct>W2 )  yield  different  MDr  .  For  example,  the  MDR 
values  based  on  the  dCo/dM  =  0.1  method  for  Cd^2  at  6.8  X  106  is  Mdr  =  0.696.  The  corresponding 
Mdr  value  for  Cpw  is  Mdr  =  0.708  for  the  same  Rc. 

Figures  17a  and  17b  show  a  similar  comparison  of  CdW2  and  Cdw  for  CL  =  0.5.  Once  again 
there  is  a  similar  behaviour.  CpW2  again  shows  a  mild  peak  at  0.011  for  Rc  =  6.8  X  106  and  =  0.64 
followed  by  a  moderate  bucket  with  CdW2  minimum  of  0.0104  at  =  0.662.  Even  at  the  higher 
Reynolds  number  Rc  =  12.8  X  106  there  is  a  vague  gesture  of  a  bucket  with  CdW2  minimum  of  about 
0.0114  at  =  °-68.  The  versus  Moo  graph  at  CL  =  0.5  in  Figure  17b  shows  a  gentle  decrease 
up  to  0.6  (up  to  0.66  for  Rc  =  6.8  X  106).  All  three  Reynolds  numbers  tested  confirm  the  bucket 
type  decrease  in  drag.  Cqw  for  the  Mach  number  range  0.66  <  Moo  ^  0.72.  A  minimum  value  of 
CDw  =  0.0057  was  observed  at  M^  =  0.68  for  Rc  =  6.8  X  106.  Once  again  the  drag  rise  Mach  number 
values  from  both  methods  are  shown  in  the  figures. 

This  theme  of  two  types  of  drag  (CDw2  and  CDw)  comparison  is  continued  in  Figures  18a 
and  18b  respectively  for  a  CL  =  0.6.  CDui  in  Figure  18a,  after  an  upward  trend  against  Mach 
number  is  showing  more  of  a  pronounced  bucket  for  Rc  =  6.8  X  106  in  the  Mach  number  range 
0.64  <  Moo  <  0.71  with  a  minimum  CdW2  value  of  about  0.0109.  The  higher  Reynolds  number, 

Rc  =  12.8  X  106,  curve  however  shows  somewhat  of  an  inflexion  region  at  Moo  =  0.66  and 
CdW2  =  0.0130  and  continues  to  grow  for  higher  Mach  number  values.  Also  included  with  the  curves 
of  Figure  18a  are  the  results  from  the  BGK  computer  code  for  Rc  =  15  X  106  which  do  not  show  any 
bucketing  which  is  to  be  expected  since  transition  is  fixed  at  7  and  15%  on  the  upper  and  lower 
surfaces.  The  corresponding  drag  results  Cpw  are  shown  in  Figure  18b.  There  is  little  change  in  drag 
up  to  Mqo  =  0.6.  The  higher  Reynolds  number  Rc  =  12.8  X  106  curve  shows  some  peaking  at  M^  =  0.6 
before  dropping  with  the  ‘bucket’  minimum  value  of  0.0082  at  M^  =  0.68.  But  the  low  Reynolds 
number  curve  continues  to  diminish  beyond  M^  =  0.6  with  the  rate  of  decrease  intensifying  after 
Mqo  =  as  it;  enters  the  bucket  behaviour  to  give  a  Cdw  =  0.0064  at  Mqo  =  0.695.  This  figure  also 
shows  differences  between  CDw  and  CdW3  which  are  typical  also  at  other  conditions. 

Note  that  the  difference  in  drag  rise  Mach  number  Mdr,  as  given  by  the  two  methods  is 
very  small.  It  lies  in  the  limits  0.7  <  Mdr  <  0.704  based  on  9Cd/9M  =  0.1  method  and  in  the  range 
0.706  <  Mdr  <  0.71  based  on  20  counts  approach. 

Figures  19a  and  19b  compare  CdW2  and  CDw  at  a  CL  of  0.7.  The  overall  drag  behaviour  in 
different  Mach  number  regimes  is  very  similar  to  the  previous  cases.  Note  the  uptrend  in  Cdw  at 
Mqq  =  0.64  before  entering  the  bucket  in  the  Mach  boundary  0.64  <  Mm  <  0.7.  A  minimum  Cdw 
value  of  about  0.0072  was  recorded  at  the  bottom  of  the  bucket  for  Rc  =  6.8  X  106  at  M^  =  0.68. 

At  Rc  =  12.8  X  106,  this  value  of  minimum  CUw  has  risen  to  0.0085  for  the  same  Mach  number 
Moo  =  0.68. 

For  CdW2,  Figure  19a  shows  the  presence  of  a  drag  bucket  in  the  region  0.64  <  M^  <  0.7, 
after  showing  a  peak  at  0.64.  The  minimum  value  of  drag  in  the  bucket  is  about  0.0114  at 
Mqq  =  0.664.  Cdw?  does  not  show  any  buckets  at  higher  Reynold  numbers. 

An  explanation  of  the  drag  buckets  mentioned  above  will  be  made  when  we  study  pressure 
distributions  in  a  later  section. 

Finally  we  show  values  of  the  range  parameter  M^Ci/Coy  on  Figure  20.  It  is  interesting 
to  note  that  for  all  three  Reynolds  numbers  the  maximum  range  parameter  is  predicted  at  a  higher  Cl 
than  design.  The  maximum  is  at  M,*,  =  0.68  and  0.7  <  Cl  <  0.8  approximately. 

Note  that  surface  roughness  does  not  seem  to  have  been  a  factor  in  triggering  transition  as 
our  surface  roughness  has  an  RMS  value  of  about  15  X  10"6  inches  compared  to  a  laminar  boundary 
layer  displacement  thickness  of  about  lO""*  inches  at  mid  chord. 
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Table  1  summarizes  the  findings  of  this  section.  Note  that  the  drag  values  are  exceptionally 
good  compared  to  other  airfoil  sections  particularly  in  view  of  the  fact  that  this  foil  is  21%  maximum 
thickness.  A  comparison  with  other  airfoils  will  now  be  made. 

4.2  Drag  Comparison  Against  Other  Airfoils 

The  drag  of  the  current  21%  thickness  airfoil  is  significantly  lower  than  other  foils  tested 
at  NAE  (except  for  the  16%  thickness  foil  from  the  same  family  of  foils).  On  inspection  of  Figure  21a 
we  see  that  Hoemer’s  (Ref.  16)  fully  turbulent  shock  free  drag  curve  provides  a  representative  low 
boundary  for  most  of  the  foils.  Note  that  some  of  the  foils  used  fixed  transition  strips  and  thus  lost 
the  advantage  of  natural  laminar  flow  (NLF).  However,  some  of  the  foils  using  NLF  still  showed  drag 
levels  comparable  with  those  of  fixed  transition.  This  could  have  been  due  to  unfavourable  pressure 
gradients  (as  in  the  ‘peaky’  type  airfoils)  or  could  be  due  to  tunnel  turbulence  levels  at  the  time  of 
testing.  Recent  improvements  (Ref.  17)  to  the  NAE  5  ft  X  5  ft  wind  tunnel  might  have  cleaned  up 
some  of  this  turbulence.  Thus  our  current  16%  and  21%  foils  may  be  taking  fuller  advantage  of  NLF. 

Also  shown  on  Figures  21a  and  21b  are  18%  and  24%  t/c  augmentor  wing  (multi-element 
airfoil)  data  as  taken  from  Reference  18. 

It  can  be  seen  that  both  16%  and  21%  foils  show  excellent  drag  levels  at  R  =  8  X  106/ft 
while  at  15  X  106/ft  the  21%  still  performs  remarkably  well  and  the  16%  is  very  good.  Even  at 
20  X  106/ft  the  21%  is  only  just  giving  a  drag  value  above  that  of  Hoemer’s. 

Other  comparisons  of  the  current  21%  airfoil  (as  well  as  our  16%)  are  made  in  Reference  1. 
In  that  reference  it  is  demonstrated  that  our  airfoil,  at  design  Mach  number,  has  drag  levels  comparable 
to  NASA  NLF  airfoils  at  Mm  =  0.1. 

Besides  comparing  the  NLF  drag  levels  with  other  airfoils  it  is  also  useful  to  compare  the 
turbulent  drag  levels  recorded  on  probe  2,  namely  Cd^-  This  value  may  be  more  representative  of 
true  flight  conditions  and  so  ideally  we  require  a  low  drag  level  here  also.  On  inspecting  Figure  21b 
we  see  that  our  drag  levels  are  very  good  compared  to  other  single  and  multi-element  airfoils  with 
the  8  and  20  X  106  Reynolds  number  values  only  just  above  Hoemer’s  curve. 


5.0  PRESSURE  DISTRIBUTIONS 

5.1  A  Comparison  With  Theoretical  Computations 

Comparison  of  experimental  and  theoretical  pressure  distributions  can  only  be  meaningful 
in  our  case  if  we  fix  transition  in  the  calculation  quite  near  to  the  leading  edge.  This  has  to  be  done 
since,  as  mentioned,  the  pressure  holdes  themselves  cause  turbulence  and  lie  in  a  turbulent  strip  of 
the  otherwise  laminar  flow  airfoil.  The  theory  used  here  is  the  BGK  non-conservative  code*9'  with 
Green’s  boundary  layer  method.  A  comparison  near  the  design  point  is  made  in  Figure  22  and  shows 
a  reasonably  good  agreement.  The  next  figure l 23 1  shows  the  effect  of  moving  the  transition  point  on 
both  the  upper  and  lower  surfaces  back  to  near  the  minimum  pressure  values  on  each  surface.  It  can 
be  seen  that  the  pressure  distributions  are  quite  different  and  that  the  pressure  becomes  more  favour¬ 
able  as  the  transition  point  moves  back  while  at  the  same  time  i  ie  shock  strength  remains  constant 
with  Moofshock)  =  1.17.  The  aft  loading  is  increased.  This  increase  in  aft  loading  will  be  investigated 
in  future  tests  when  some  pressure  holes  aft  of  60%  will  be  placed  at  a  different  spanwise  location  on 
the  airfoil. 


The  drag  values,  for  =  0.68  and  Cl  =  0.6,  for  different  theoretical  locations  of  transition 
are  shown  in  Figures  24a  and  24b.  As  expected  the  drag  differences  are  very  significant.  Our  experi¬ 
mental  drag  is  also  shown  as  a  locus  on  the  figures.  As  can  be  seen,  correlation  of  experimental  drag 
with  the  theoretical  values  obtained  with  natural  transition  is  difficult  as  the  higher  Reynolds  number 
case  yields  a  lower  drag  (0.0045)  than  the  lower  Reynolds  number  case  (0.0057).  TTiis  is  opposite  to 
the  experimental  observation  (0.0082  and  0.0071). 
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Another  attempt  at  a  correlation  can  be  made  if  we  adjust  upward  the  theoretical  drag 
levels  by  22  counts  (found  from  previous  correlations  on  airfoils  with  smaller  runs  of  laminar  flow). 

To  observe  this  correlation  we  have  plotted  on  Figures  24  the  locus  of  the  experimental  drag  levels 
minus  22  counts.  It  can  be  seen  from  Figure  24b,  at  Rc  =  6.8  X  106,  that  a  reasonable  prediction 
of  drag  could  be  made  by  assuming  about  65%  and  40%  transition  points  for  the  upper  and  lower 
surfaces  respectively.  At  Rc  =  12.8  X  106,  Figure  24a,  the  indication  is  that  transition  would  be  closer 
to  50%  and  30%  respectively. 

5.2  Experimental  Pressure  Behaviour 

The  effect  of  Reynolds  number  on  the  pressure  distribution  is  shown  in  Figure  25. 
Remembering  that  these  pressure  taps  are  in  a  turbulent  boundary  layer  it  is  not  surprising  that  the 
pressures  are  virtually  independent  of  Reynolds  number.  One  can  only  surmise  that  had  the  pressure 
measurements  been  taken  in  the  presence  of  the  natural  laminar  boundary  layer  the  distributions 
would  have  been  different  due  to  the  different  lengths  of  laminar  flow. 

The  difference  in  pressure  distributions  as  Mach  number  is  increased  for  a  constant  lift 
of  roughly  0.65  is  shown  in  Figure  26.  Assuming  that  at  least  the  trends  of  this  turbulent  pressure 
distribution  are  similar  to  the  NLF  trends  we  can  conclude  that  as  increases  the  pressure  is 
becoming  more  favourable  to  longer  runs  of  laminar  flow.  This  must  be  the  reason  for  the  drag  bucket 
mentioned  earlier. 


6.0  CONCLUSIONS 

At  the  lower  Reynolds  number  of  6.8  X  106  it  appears  that  long  runs  of  laminar  flow  are 
possible  producing  very  low  drag  values.  The  length  of  this  flow  diminishes  as  Reynolds  number  is 
increased  until  at  16.7  X  106  drag  levels  are  comparable  to  other  “turbulent  flow”  airfoils  tested  at 
NAE. 


Because  of  the  success  of  this  airfoil  and  of  a  similar  16%  design  (Ref.  1)  further 
investigations  will  be  made  on  thinner  airfoil  sections  which  will  exploit  natural  laminar  flow. 
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TABLE 1 

A  SUMMARY  OF  PARAMETERS  CDw  AND  CDwz 
CDw 


CONDITIONS 

10“®  Rc 

CL  =  0.3 

■531 

CL  =  0.6 

CL  =  0.7 

Bucket  Region 

7 

0.66  <  M^  0.712 

0.66  <M„,<  0.712 

0.66  <M_t<  0.710 

0.64  <Moo<  0.710 

13 

0.66  <Moo<  0.718 

0.66  <M^  <0.710 

0.66  <M  <0.700 

- 

17 

0.66  <  M^  0.700 

0.66  0.696 

— 

Bucket  Min  Drag 

7 

0.0057 

0.0057 

0.0064 

0.0072 

13 

0.0065 

0.0079 

0.0082 

0.0085 

17 

0.0099 

0.0099 

0.0110 

— 

Mjm  based  on  aC^/aM^  =  0.1 

7 

0.708 

0.703 

0.700 

0.692 

13 

0.704 

0.703 

0.687 

0.688 

17 

— 

0.699 

— 

0.690 

Mdr  based  on  20  counts 

7 

0.720 

0.719 

0.710 

— 

13 

0.722 

0.716 

0.700 

0.700 

17 

— 

~~ 

— 

— 

Cl>W2 


CONDITIONS 

ur*  rc 

CL  =  0.5 

CL  =  0.6 

Ch  =  0.7 

Bucket  Region 

7 

0.64  <Moo<  0.700 

0.64  <Mw<  0.686 

0.64  <  Moo<  0.704 

0.64  <M„<  0.696 

13 

17 

_ 

0.66  <  Ms  <0.686 

_ 

Bucket  Min  Drag 

7 

0.0110 

0.0104 

0.0109 

0.0114 

13 

17 

— 

0.0114 

_ 

— 

Mob  based  on  dCu/aM^  =  0.1 

7 

0.696 

0.694 

0.704 

0.699 

13 

17 

0.662 

0.690 

0.700 

0.690 

Mqr  based  on  20  counts 

7 

0.702 

0.694 

0.706 

13 

0.624 

0.690 

0.657 

0.682 

17 

— 

— 

— 

— 

•  «  '' 

•  *  -  •  »  ■  .  ■ 
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FIG.  5:  CL  AND  CM  VERSUS  a,  Mm  =  0.3,  Rc  ^  6.8  AND  12.8  X  106 


FIG.  8:  LIFT  PERFORMANCE  AT  R„  =  6.8  X  106 


FIG.  10:  THE  WAKE  RAKE  PROBE  LOCATIONS  RELATIVE  TO  THE  AIRFOIL 
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SW  RUN 

o 

CYC  M 

Rcxl0“6 

V/U 

°  28369 

/  1 

1  0.801 

12.76 

0.0084 

o  28377 

/  1 

1  0.600 

6.78 

0.0084 

*  28377 

/  2 

1  0.602 

6.87 

0.0083 

BALANCE 

DATA 

WITH  WAKE 

DRAG 

PASS  2  DATA  - 

CORRECTED 

SYM  RUN 

CYC  A 

Rcxl0'6 

V/U 

°  28388 

/  1 

1  0.700 

12.72 

0.0084 

=  28388 

/  2 

1  0.700 

12.82 

0.0083 

*  28380 

/  1 

I  0.700 

6.79 

0.008] 

•  28380 

/  2 

1  0.701 

6.85 

0.0080 

BALANCE 

DATA 

WITH  WAKE 

DRAG 

{  PASS  2  DATA 

-  CORRECTED 

0.028 


0.024 


0.020 

CD 

0.*0!6 


0.012 


0.008 


0.004 


FIG.  13a:  CDvy  VERSUS  CLb  FOR  M 
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SYM  RUN  CYC  M 

Rcxl0~6 

V/U 

»  28367  / 1 

!  0.301 

12.69 

0.0081 

o  28368  /] 

1  0.498 

12.69 

0.0084 

*  28384  / 1 

1  0.600 

12.69 

0.0084 

•  28384  /2 

1  0.600 

12.83 

0.0083 

BALANCE  DATA 

WITH  WAKE  DRAG 

PASS  2  DATA  - 

CORRECTED 

3$: 


-  0.004 


III! 


FIG.  14b:  CDw  VERSUS  CL_  FOR  =  0.3,  0.5  AND  0.6 


SYM  RUN 

CYC  M 

Rcx  1 0 

6  V/U 

°  28385 

/  1 

1  0.680 

12.76 

0.0083 

o  28385 

/  2 

1  0.681 

12.87 

0.0082 

*  28388 

/  1 

1  0.700 

12.72 

0.0084 

•  28388 

/  2 

1  0.700 

12.82 

0.0083 

«  28389 

/  1 

1  0.661 

12.69 

0.0083 

»  28389 

/  2 

1  0.661 

12.78 

0.0083 

BALANCE 

DATA 

WITH  WAKE 

DRAG 

PASS  2  DATA 

-  CORRECTED 

VERSUS  M„  FOR  C 


MACH  NUMBER  SHOWN  AS  INTEGER 

ALL  AIRFOILS  TESTED  AT  NAE  5  FT  x  5  FT  FACILITY 


‘G  COMPARISONS  WITH  OTHER  AIRFOILS  USING  C, 


MACH  NUMBER  SHOWN  AS  INTEGER 

ALL  AIRFOILS  TESTED  AT  NAE  5  FT  X  5  FT  FACILITY 


FIG.  23:  TO  SHOW  THE  EFFECT  OF  DIFFERENT  TRANSITION  POINTS  ON  THE 
BGK  PRESSURE  DISTRIBUTION  Mm  =  0.68,  CL  =  0.623,  Rc  =  12.8  X  106 


LOCUS  OF  C0=  0.0082 


LOWER 
SURFACE 
TRANSITION 
(%  x/c) 


84^^  •  74 

% 

% 


•91  \  *70  > 

\  67* 
86»\ 

\  *79 


(NATURAL  TRANSITION- COMPUTED) 


-LOCUS  OF  CD-  0.0022 

w 


>75  *65 


UPPER  SURFACE  TRANSITION  (%x/c) 

FIG.  24a:  DRAG  COUNT  LEVELS  FOR  VARIOUS  POSITIONS  OF  TRANSITION 
(THEORETICAL)  Rc  =  12.8  X  106,  =  0.68,  CL  =  0.6 

EXPERIMENTAL  Cqw  =  0.0082 


(NATURAL  TRANSITION- 
COMPUTED)  \ 


LOWER 

SURFACE  An 
TRANSITION 
(%  x/c) 


9  *87  ^L*68 

*96 

)6  *83  V, 


-LOCUS  OF  CD-  0.0022 

w 


LOCUS  OF  CD=  0.0071 


•87  • 77 


UPPER  SURFACE  TRANSITION  (%  x/c) 


FIG.  24b:  DRAG  COUNT  LEVELS  FOR  VARIOUS  POSITIONS  OF  TRANSITION 
(THEORETICAL)  Rc  =  6.8  X  106,  M*,  =  0.68,  Cu  =  0.6 
EXPERIMENTAL  Cr»...  -  0.0071 


EASING  Mm  ON  THE  Cp  DISTRIBUTION 
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